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Abstract 

By applying the closed-time-path Green function formalism to the chiral dynamical model based 
on an effective Lagrangian of chiral quarks with the nonlinear-realized meson fields as bosonized 
auxiliary fields, we then arrive at a chiral thermodynamic model for the meson fields after integrat- 
ing out the quark fields. Particular attention is paid to the spontaneous chiral symmetry breaking 
and restoration from the dynamically generated effective composite Higgs potential of meson fields 
0^ at finite temperature. It is shown that the minimal condition of the effective composite Higgs 

p) I potential of meson fields leads to the thermodynamic gap equation at finite temperature, which 

^ ■ enables us to investigate the critical behavior of the effective chiral thermodynamical model and 

to explore the QCD phase transition. After fixing the free parameters in the effective chiral La- 
grangian at low energies with zero temperature, we determine the critical temperature of the chiral 
symmetry restoration and present a consistent prediction for the thermodynamical behavior of 
several physically interesting quantities, which include the vacuum expectation value Vo{T), quark 
condensate < qq > (T), pion decay constant f^iT) and pion meson mass mT^{T). In particular, it 
is shown that the thermodynamic scaling behavior of these quantities becomes the same near the 
critical point of phase transition. 
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I. INTRODUCTION 



Thermodynamics of quantum chromodynamics (QCD) has been attracted a lot of at- 
tention during the last three decades. Many interesting physical phenomena are related to 
it, such as the equation of state of quark gluon plasma(QGP), chiral symmetry breaking 
and restoration, deconfinement phase transition, and so on. Thus, the study of the QCD 
thermodynamics becomes a basic problem in our understanding of the strong interaction. 
Particularly, the deconfined QGP is expected to be formed in ultrarelativistic heavy-ion col- 
lisions(HIC) [ll-Isj and many experiments like those at the Relativistic Heavy Ion Collision 
(RHIC) and at the Large Hadron Collider (LHC) have been built to explore the nature of 
the QGP and to search for the critical point of the phase transition, which enables many 
theoretical ideas testable and makes the research area more exciting. 

In this paper, we are going to study the thermodynamic properties of the effective chiral 
dynamical model(CDM)[9j of low energy QCD with spontaneous symmetry breaking via 
the dynamically generated effective composite Higgs potential of meson fields. The main 
assumption in such an effective CDM was based on an effective Lagrangian of chiral quarks 
with effective nonlinear-realized meson fields as the bosonized auxiliary fields, which may 
be resulted from the Nambu-Jona-Lasinio(NJL) four quark interaction 10| due to the strong 
interactions of nonperturbative QCD at the low energy scale, then after considering the 
quantum loop contributions of quark fields by integrating over the quark fields with the loop 
regularization (LORE) method|lll.ll2|. the resulting effective chiral Lagrangian for the meson 
fields in the CDM@ has been turned out to provide a dynamically generated spontaneous 
symmetry breaking mechanism for the SU{3)l x SU{3)Fi chiral symmetry. The key point 
for deriving the dynamically generated spontaneous symmetry breaking mechanism is the 
use of the LORE method which keeps the physically meaningful finite quadratic term and 
meanwhile preserves the symmetries of original theory. More specifically, the advantage of 
the LORE method is the introduction of two intrinsic energy scales without spoiling the basic 
symmetries of original theory, such intrinsic energy scales play the role of the characteristic 
energy scale Mc and the sliding energy scale fig- Here Mc is the characterizing energy scale of 
nonperturbative QCD below which the effective quantum field theory becomes meaningful 
to describe the low energy dynamics of QCD, and /i^ reflects the low energy scale of QCD 
on which the interesting physics processes are concerned. As a consequence, it was shown 
that the resulting effective CDM of low energy QCD can lead to the consistent predictions 
for the light quark masses, quark condensate, pseudoscalar meson masses and the lowest 
nonet scalar meson masses as well as their mixing at the leading order [3]- Based on the 
success of the effective CDM at zero temperature, we will show in this paper how the CDM 
can be extended to an effective chiral thermodynamic model (CTDM) at finite temperature 
by applying the closed-time-path Green function(CTPGF) approach (l4 - 19| . and how the 
CTDM enables us to describe the critical behavior of the low energy dynamics of QCD, 
and in particular to determine the critical temperature of the chiral symmetry breaking 
and restoration. This comes to our main motivation in the present paper. For our present 
purpose, we will investigate the CTDM with considering only two flavors and ignoring the 
possible instanton effect which is known to account for the the anomalous U{1)a symmetry 
breaking. 

Following the almost same procedure in deriving the CDM for the composite meson fields, 
we will arrive at an effective CTDM at finite temperature by adopting the CTPGF formalism. 
The key step in the derivation is the replacement of the zero-temperature propagator for the 
quark field with its finite temperature counterpart in the CTPGF formalism which has been 
shown to be more suitable for characterizing the nonequilibrium statistical pro cesses [l6k 



19j. We then obtain the dynamically generated effective composite Higgs potential at finite 
temperature, its minimal condition leads to the gap equation at finite temperature and 
enables us to explore the critical behavior and temperature for the chiral symmetry breaking 
and restoration. After fixing the free parameters in the effective Lagrangian, we present our 
numerical predictions for the temperature dependence of physically interesting quantities, 
such as the vacuum expectation value Vo(T), the pion decay constant fniT), and the masses 
of pseudoscalar and scalar mesons. The resulting critical temperature is found to be Tc ^ 200 
MeV, which is consistent with the NJL model prediction |20l-l23| . 



II. OUTLINE ON CHIRAL DYNAMICAL MODEL OF LOW ENERGY QCD 

Before exploring the thermodjTiamic behavior of the chiral dynamical model, it is useful 
to have a brief review on its derivation of the effective chiral Lagrangian at zero-temperature. 
For our present purpose with paying attention to the investigation of the chiral symmetry 
restoration at finite temperature, we shall consider the simple case with only two flavors 
of light quarks u and d, and ignore the instanton effects. For simplicity, we also assume 
the exact isospin symmetry for two light quarks q = {u, d)'^ by taking the equal mass 
rriu = rrid = m. Our discussion here is mainly following the previous paper by Dai and 
Wu|9|]. 

Let us begin with the QCD Langrangian for two light quarks 

Cqcd = qi'{id, + gsG1T'^)q - qMq - UrG.^G^'' (1) 

where q = {u, d)"^ denotes the SU(2) doublet of two light quarks and the summation 
over color degrees of freedom is understood. are the gluon fields with SU(3) gauge 
symmetry and Qs is the running coupling constant. M is the light quark mass matrix 
M = diag{mi,m2) = diag{mu-,md). In the limit mj — )■ (i = 1,2), the Lagrangian has the 
global U{2)l X U{2)^ chiral symmetry. Due to the instanton effect, the chiral symmetry 
U{1)l X U{1)r is broken down to the diagonal U{l)v symmetry. This instanton effect can 



be expressed by the effective interaction [24, |25| 

C""'' = t^inste''-' deti-qmi) + h.c. (2) 

where t^inst is the constant containing the factor e"^'^^/^^. Obviously, such an instanton term 
breaks the U{1)a chiral symmetry. As our present consideration is paid to the phenomenon 
of the SU{2)l X SU{2)r chiral symmetry breaking and its restoration at finite temperature, 
we will switch off the instanton effect by simply setting Hmst = in the following discussions. 

The basic assumption of the CDM is that at the chiral symmetry breaking scale (~ IGeV) 
the effective Lagrangian contains not only the quark fields but also the effective meson fields 
describing bound states of strong interactions of gluons and quarks. After integrating over 
the gluon fields at high energy scales, the effective Lagrangian at low energy scale is expected 
to have the following general form when keeping only the lowest order nontrivial terms: 

£e//(g, g, = qi^id^^q + qLl^,A1qL + qRli^A^RqR - [qd'^ - M)qR + h.c] 

+nl^tr{<l>M^ + M$"^) - /ijtr$$"^ (3) 

where $jj are the effective meson fields which basically correspond to the composite operators 
qRjqhi- Al and Ar are introduced as the external source fields. It is noticed in Eq.Q that 



the effective meson fields $jj are the auxihary fields in the sense that there is no kinetic term 
for them, which may explicitly be seen by integrating out we then obtain the following 
effective Lagrangian of quarks: 



- l){qLMqR + quM^qL) + —qLqmRqL (4) 



which arrives exactly at the Nambu-Jona-Lasinio(NJL) model[lO| of effective four-quark 
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interaction with the quark mass matrix — 1)M. When the mass term is understood as 

the well-defined current quark mass term in the QCD Lagrangian Eq.([T]), it is then clear 
that /i^/ = 2. Note that the high order terms with the dimension above two for the 
effective meson fields are not included, which are assumed to be small and generated in 
loop diagrams, so only the lowest order nontrivial fermionic interaction terms are taken into 
account after integrating over the gluon field. 

For our present purpose, we only focus on the scalar and pseudoscalar mesons while the 
vector and axial vector sectors will not be considered here. As the pseudoscalar mesons are 
known to be the would-be Goldstone bosons, the effective chiral field theory is naturally to 
be realized as a nonlinear model. Thus we may express the effective meson fields $(a;) into 
the following 2x2 complex matrix form: 

$(x) = U{x) ^ = = 

3 3 

(t)\x) = (f){x) = J2 (t>''{x)T\ I{\x) = U{x) = n"(x)T" (5) 

a=0 a=0 

where (a = 0, 1, 2, 3) with [T", T*] = i/'^'"^T^ and 2trT''T^ = dab are the four generators of 
f/(2) group. The fields 11'* (x) represent the pseudoscalar mesons and 0°(x) the corresponding 
scalar chiral partners. Where / is known as the decay constant with mass dimension. 

The Lagrangian Eq.dH]) with the definition of composite meson fields Eq.dS]) is 

our starting point for the derivation of the effective chiral Lagrangian for mesons. The 
Lagrangian can be obtained by integrating over the quark fields and the procedure for the 
derivation can be formally expressed in terms of the generating functionals via the following 
relations 



Let us first demonstrate the derivation of Ceff{^). In order to obtain the effective chiral 
Lagrangian for the meson fields, we need to integrate over the quark fields (which is equiva- 
lent to calculate the Feynman diagrams of quark loops) from the following chiral Lagrangian 

^Iffil, q) = qi^id^q + qili^A^qL + qRlt^A^RqR - [qd^ - M)qR + h.c] + xm (7) 

where we have introduced the real source field x(^) for the composite operator qq and the 
source term XQQ which will be useful for the derivation of the chiral thermodynamic model, 
while eventually the source field is taken to be zero x = 0. 



With the method of path integral, the effective Lagrangian of the meson fields is evaluated 
by integrating over the quark fields 

[d<^>]exp{i J d^xC^] = Z^^ J [dmdq][dq]exp{i J d^xCl^^} (8) 

The functional integral of the right hand side is known as the determination of the Dirac 
operator 

J [dq][dq]exp{i J d^xCl^^} = det(iP^) (9) 
To obtain the effective action, it is useful to go to Euclidean space via the Wick rotation 

7o ilA-i Go iG^-i — >■ —1x4^ (10) 
and to define the Hermitian operator 

Se ^ j d^XEC^ = In det{iVl) 

= ^[lndet(iP|) + lndet(iP|)^] + ^[Indet(iP^) - lndet(iP|)^] 

with 

S^ne = I (^^ECg = ^ In det(iP|(iP|)t) ^ 1 In det A| - In Zo (12) 

S^im = / d'xECr^ = ^ In det(iP|/(iP|)t) ^ 1 In det (13) 

where the imaginary part £^ appears as a phase which is related to the anomalous terms 
and will not be discussed in the present paper. The operators in the Euclidean space are 
given by 

= -i7 ■ a - 7 • AlPl - 7 • ArPr + ^Pr + $^Pl + X 

= i^E + X 

(iP|)t = i7 . 9 + 7 ■ ArPl + 7 • AlPr + ^^Pr + ^Pl + X 

= {iVe)^ + X (14) 

with $ = $ - M and P± = (1 ± 75)72. iV\, {iV^Y and A| are regarded as matrices in 
coordinate space, internal symmetry space and spin space. Noticing the following identity 

IndetO^TrlnO (15) 

with Tr being understood as the trace defined via 

TrO ^tr d^x < x\0\y > \a,=y (16) 



Here tr is the trace for the internal symmetry space and < a;|0||/ > is the coordinate matrix 
element defined as 

(17) 



For the derivative operator, one has in the coordinate space 

^x\d^'\y>=6\x-y){-ik^' + d^) (18) 

With the above definitions, the operators i'D'^, {iTf^Y and A'^ in the Euchdean space are 
given by 

zV'k = -^-k-i^-d-^- AlPl - 7 • ArPr + ^Pr + ^^Pl + X 

= --i-k + iV\ (19) 

ilV),)^ = ^.k + l^-d + ^-ARPL + l-ALPR + ^^PR + ^Pl + X 

= 7 ■ + {iVl)^ (20) 
= k'^ + WPR + ^^^PL-i-i ■ De^Pl-^I ■ De^^Pr 

-CTf^uJ^Rf^uPL - CTf^uJ^Lf^uPR + {iD Et,){iD Et,) + 2k ■ {iDe) 

■ dx + x\iT^E + {iVe)^] + 
= e + ^1 = ^ + Ae + A^ (21) 



where 



iDe^ = id<^ + - (22) 
iDe = td + ArPl + APij (23) 



and 



Ae = M^Pr + ^^^Pl-i-t ■ De^Pl-ii ■ De^^Pr 

-(r^,uJ^Rf,^PL - (yt,uJ^Lp,yPR + iiDE^)iiDE^) + 2k ■ {IDe) (24) 

= -1-1 ■ dx + xVT^E + {^e)^] + (25) 

A^ = As + A^ (26) 

Thus the effective action is obtained as 
NT r fl^k 

where the subscripts SF refer to the trace over the spin and fiavor indices and is the 
color number. 

Before proceeding, we would like to mention some formulae which will be useful for the 
derivation of chiral thermodynamic model late on. By taking the functional derivative of 
Eq.f l27|) with respect to the source field x at x = 0, we have 



^ J d^XEtrsF J (-7 ■ ^ + 7 ■ (^ii - A)75 + $ + 



5xix) 



\x=o 



k^ + AE 



from the right hand side one may pick up the quark propagator as x is the source field 
for the quark operator q{x)q{x). This can explicitly be shown by differentiating the ef- 
fective action with respect to x? which gives the coinciding limit of the quark propagator 
lim^^^j^ trsF{T[qix)q{y)]). 



Alternatively, if taking the functional derivative of Eq.( l27|) with respect to the operator 
A^, we obtain another form of the formulae 



SiA^)ij 2 J (27r)4 yk^ + AE + A 



(29) 



X/ ji 



In other word, when functionally integrating over A-^ and summing over the flavor and spin 
degrees of freedom, we are led to the effective action Eg. ( 1271) . The physical action is yielded 
by taking the source field to be zero x(a^) = 

To derive the effective chiral Lagrangian for the meson fields, we shall make the following 
redefinition for A^ = k"^ + Ae 

A| = Ao + Ae (31) 

Ao = P + M^ (32) 
Ae = ($1'"^ - M^)Pr + (<l"f<l - M^)Pl - ^7 ■ De^Pl - il ■ De^^Pr 

-(^f^uJ^Rf^uPi - (yfiu^Lf^uPR + iiDEp,)iiDEp,) + 2k ■ {iDe) (33) 

where M is the supposed vacuum expectation values (VEVs) of $, i.e., < $ >= M = 
diag.{ffiuifh(i) with ffiu = fhd = fh under the exact isospin symmetry. Here rhi = Vi — rrii is 
regarded as the dynamical quark masses, and Vi is supposed to be the VEVs of the scalar 
fields, i.e., < >= V = diag. {vi,V2) which will be determined from the minimal conditions 
of the effective potential in the effective chiral Lagrangian Ceff{^)- With this convention, 
it is seen that the minimal conditions of the effective potential are completely determined 
by the lowest order terms up to the dimension four ($$"'" — M^)^ in the effective chiral field 
theory of mesons. _ 

By regarding A^; as the perturbative interaction term and taking Zq = (det Ao)2, the 
effective action in the Euclidean space can be written as 

AT r r d^k 

S^Re = Yj ^'^^ J -^,trsF[H^o + AE)-\nAo] 
Nc f ,A f d^k , 1 ~ , 

2 J (2ir)4 ^'^^ n ^ A„ 

[ ,. [ d*k , 1 J 11 J,, 

Note that we only keep the first two terms in the expansion over A^; since these are the 
only divergent terms in the integration over the internal momentum k. Particularly, the 
divergence degree of the first integral is quadratical while the second logarithmic. In order 
to maintain the gauge invariance and meanwhile keep the divergence behavior of the integral. 



we adopt the loop regularization (LORE) method proposed in (ill . |12| for the momentum 



integral 



d^k 1 M2 /i 



2 



(27r)4 F + M2 ^ ^2 - (35) 
' ^^^ = T^^o(^) (36) 



(27r)4 (P + M2)2 " 167r2 "^M2 
with the consistent conditions for the tensor type divergent integrals 



^2niy ~ 2^^^'^^^'' "^oa*!^ ~ ^9^iuIq (37) 
where 



d k kf_i ky J C d J\i 



' (27r)4(P + M2)2' 0'^'^ J (27r)4 (F + M2)3 



(38) 



with 



The two diagonal matrices Lq = diag.{L^\ L^^) and L2 = diag.{L^\ L^^) are given by the 
following form: 

4..^,„M_,„^,„,^) (39) 

/•^ 1 _ g-'^ 1 

Voix) = / dcr , yi{x) = -{e ^ -1 + x), y2{x) = yo{x) - yi{x) (41) 

Jo ^ 

Note that is the characteristic energy scale from which the nonperturbative QCD effects 
start to play an important role and the effective chiral field theory is considered to be valid 
below the scale Mc- We have also introduced the definitions 

f^l = f4 + "^i; Thi = Vi- rrii (42) 

with fil the sliding energy scale. It is usually taken to be at the energy scale at which the 
physical processes take place, which is expected to be around the QCD scale Aqcd for our 
present consideration. 

With these analysis, the effective chiral Lagrangian can be systematically obtained to be 



S^Re = ^/ d'xEtMM'^L2[m^-M') + {^^^-M')] 



~Lo[De^ ■ De^^ + De^^ ■ De^ + ($$"^ - M^f + (<|)"f<l - M^f]} (43) 

where the trace over the spin indices give the factor 2 since our quark fields defined in Eq. ([3]) 
are Weyl fermion fields, each of which has 2 degrees of freedom. 

By transforming back to the Minkowski spacetime signature and adding the extra terms 
in Eq.(|3]), we finally arrive at the following effective chiral Lagrangian at zero temperature 

+ ^M>fL2[(<l>^<i> - M') + (Ixl^ - Af2)] 

+/i^trir($M"f + M$"f) - jjjtr^^l (44) 



The derivation of the above effective chiral Lagrangian in an equivalent rotated basis |1 3 
is given in Appendix |Bl which may be more transparent for the spontaneous symmetry 
breaking with the composite Higgs-hke scalar mesons. 



III. DERIVATION OF CHIRAL THERMODYNAMIC MODEL OF QCD 



After a brief outline for the derivation of the effective chiral Lagrangian of the CDM 
for mesons at zero temperature, it is straightforward to incorporate the finite temperature 
effects into the effective Lagrangian. The method for the derivation of finite temperature 
effective Lagrangian is similar by applying for the closed-time-path Green function(CTPGF) 
formalism. The CTPGF formalism, developed by Schwinger and Keldysh [Hj, has 
been used to solve lots of interesting problems in statistical physics and condensed matter 
theory It is generally believed that this technique is quite efficient in investigating the 
nonequilibrium and finite temperature dynamical systems as this formalism simultaneously 
incorporates both the statistical and dynamical properties (Tg], 
to treat a system of self-interacting bosons described by Xcf)'^ 



17 



It has also been used 
scalar fields j26|. A brief 
introduction to this formalism is given in Appendix A. Readers who are not familiar with 
the CTPGF formalism are refered to the excellent review articles fl3, 27 1 and monographs 



18, 19 



As shown in the Appendix A, the main step for the derivation of the effective action with 
finite temperature in the CTPGF formalism is to replace the field propagator with its finite 
temperature counterpart [l6|, 17, 27, 28|. Applying the CTPGF formalism to the propagator 
of the quark fields in Eq. fl28l) . we arrive at the following result 



6S 



M 



5xix] 



\x=o 



2 



1 



2m9{h)6{k^ + Ae 



2TTie{-h)6{k^ + Ae 



1 



- 2TTinF{u))S{k'^ + Ae) 



1 1 
1 1 



(45) 



where u is defined as the effective energy oj = \J k? + Ae and npiuj) represents the Fermi- 
Dirac distribution function 

1 



^Hc^) = (46) 

For our present purpose, we only need to calculate the first component of the effective action 
since it is the only one which is related to the causal propagation j27l. 28|. For convenience, we 
may adopt the following alternative formula which is similar to Eq.f l29p at zero temperature 



6A^ 



d^k 



(27r)4^P + As + A, 
d^k 1 
(27r)4 F + + 



- 2ninF{uj)S{k^ + A^ + A^)] 



d^k 



nE^uj) 



1 



k^ + AE + A, 



(47) 



By functionally integrating over A^ and summing over the spin and flavor indices, we then 
obtain the effective action for the CTDM 



d^k 



-InZo 



ln(l + e 



)] 



(48) 



where we have put the source field = in the end of the calculation. 

Separating A\ as in Eq. (l3T]) and identifying A^; as the perturbation part, we then make 
the expansion in terms of A^;. Here we only keep the lowest order terms as they are the 
only ones relevant to our present discussion. Also, we take 

In Zo = f I d^x^trs,[ j ^ ln(fc^ + + ^ / (^i ^^^^ + e'^^^)] (49) 

which provides the cancelation for the infinite zero-point energy. 

Thus, the effective action with leading terms in the Euclidean space can be written as 

r.M /■ ,4 rr /" I [ <fk 1 , ~ 

2 J J (27r)4 k^ + M^ J {2ny ^/j^^^^]^(^^i3^k^+Kp _,_ 

I, [ d^k 1 f d^k + Af2e/3Vfc2+M2 ^ ^py/P+NP + 1^X2. 

]A|} 



= / rf'a;str^{M2L2(T)[(<l>$t - M^) + ($t$ _ ^2^^ 

-^Lo(T)[L)b<|) ■ L)i5<l>"f + De^^ ■ De^ + ($$"^ - M^f + ($"^<i' - M^f]} (50) 



where we have defined 



Lo(T) = ^0(777) -- d k- 



2 



L2(T) ^ L2(^) - /" rf^fc^==^^_ (51) 

VA?2 + M2(e/3Vfc^ + 1) 



and used the results 

trsAE = 2[(<l<|t - M^) + (<|t$ _ ^2^)] (^52) 
ir^Al; = 2[L)£;<I • De^^ + /^s*!^ ■ + {^^^ - M^f + (<l"f<l - M^)^] (53) 

Finally, transforming the action to the Minkowski spacetime and adding the extra terms 
in Eq.(l3]), we arrive at the following effective chiral Lagrangian at finite temperature for the 
composite meson fields 

Ceffm = ]^-^trEL^{T)[D^^^D^'^ + D^^D^^^ -{^^^-M^f-{^^^ -M^f] 
rM>^L2(T)[(<|t<| - + ($|,t _ ]v/2)| 



167r2 

+fi'i{T)trF{^M^ + M^^) - ^J(T)tr$<l>"f (54) 

where Lq{T) and L2{T) are given in Eg. (151 p . Note that the initial mass scale fif (fim) char- 
acterizes the nonperturbative gluon effect at zero temperature. At the finite temperature, 
it is expected that the mass scale fif{fim) is in general temperature dependent, which will 
be seen more clear below. 



IV. DYNAMICAL SYMMETRY BREAKING AND THERMODYNAMIC PROP- 
ERTIES OF CTDM 



Let us now focus on the dynamically generated effective composite Higgs potential of 
meson fields, which can be reexpressed as the following general form 

+itrj.A(T)[(<l<|t)2 + (<|,t$)2] (55) 
with pif{T), fi^{T) and X{T) the three diagonal matrices 

il}{T) ^ n}{T) - ^,{MlL,{T) + M^Lo(T)) (56) 
(tiT) ^ uliT) - ^,{MlL,{T) + M^Lo(T)) (57) 

Taking the nonlinear realization $(a;) = ^L{x)(j){x)^'^j^{x) with supposing that the minimal 
of the above effective potential occurs at the point < 4> >= V{T) = diag.{vi{T),V2{T)), we 
can write the scalar fields as follows 

(P{x) = V{T) + ^{x) (59) 

where the VEVs may be written in terms of the following general form: 

Vi{T) = Voir) + (3omi i = 1,2 or i — u,d (60) 

For the equal mass m„ = md = m considered in our present case, it leads to the general 
VEVs Vi{T) = V2{T) = v{T) and the single form v{T) = VgiT) + Poin. 

By differentiating the effective composite Higgs potential of the scalar meson field at the 
VEV v{T), we then obtain the minimal conditions: 

- fx}{T)iv{Th + /i™(T),m, - 2A(T),m3(T), = (61) 

with equal mass of two flavor quarks, it reduces to one minimal condition. For convenience 
of discussions, it is helpful to decompose A*^(T), /ij(T), /i^(T) and A(T) into two parts with 
one part independent of the current quark mass m. Practically, it can be done by making 
an expansion with respect to the current quark masses 

/^^(T) = /x^(T) + 2(/3„ - l)voiT)Tfi, ^iliT) = ^l + vl{T), 

m{T) = m[l + {l3o - l)m/(2t;o(T))] (62) 

i,){T) = n}{T) + 2nfo{T)rh{T)[l + J] a,(r)(^)'=(/3, - 1)^=] (63) 

fc=i ^"'^^ 

i^liT) = flliT) + 2nf^{T)m{T)[l + '^>^(T){^^rWo - 1)1 (64) 



k=l 

x{T) = A(r) - AoX;-^'^(^)(^)'(/5o - 1)1 = 



(65) 



By keeping only the nonzero leading terms in the expansion of current quark masses, we 
then obtain the following constraints from the minimal condition Eq.f l6ip 



fij{T) + 2X{T)vi{T) = (66) 

Here the temperature-dependent parameters /i^(T), fij-iT) and A(T) are related to the initial 
parameters in the effective potential and the characteristic energy scale via the following 
relations 

/iJ(T) = /i5(T) - ^,{M'MT) + vl{T)U{T)) (67) 



-^UT) = ^^l(T) - ^(M'MT) + vliT)UiT)) (68) 

A(T) = ^MT) (69) 

where Lo{T) and L2{T) represent the leading order expansion of Lq{T) and L2{T) with 
respect to m. Explicitly, they are given by 



Lo{T) ^ Lo(^) - - / d'k-^ "-^ (70) 

UT) ^ L,(4^) - -i- / d^k ' ,^ (71) 



M2 ' 7rM2 



y^^7^l2^(e'^^^^+^ + 1) 



With these definitions of parameters, the minimal condition Eq.f l66|) is transformed into the 
following form 



Stt^ 

r,, (72) 



^0 ^A;2 + t;2(r)(e^V'='+-?(^) + l) 

which is the gap equation at finite temperature. In order to obtain the critical temperature, 
let us make a simple assumption for the temperature dependence of the mass scale fJ''f{T) 

^}{T) = ^,vl{T) (73) 

with 7o a temperature independent constant. The reason for this assumption will become 
manifest below from the thermodynamic property of the pion meson mass. In fact, recall that 
the appearance of /i^ in the effective Lagrangian Eq. ([3]) can be traced back to the integration 
over the gluon fields. In principle, it can be calculated from QCD and given in terms of 
the QCD parameters 5's(/i) and Aqqd- So the thermodynamic property of /Zj is expected 
to be obtained from the detailed analysis of gluon dynamics at finite temperature. The 
above simple assumption means that the gluon thermodynamics has the same temperature 
dependence as the chiral thermodynamics of quark condensate. 

With such an assumption and based on the chiral thermodynamic gap equation, we are 
able to calculate the critical temperature for the chiral symmetry restoration. Suppose that 



at the critical temperature the VEV Vo(T) approaches to vanish, so does the fJ.'jiT), then 
the gap equation becomes 

-7.+ l + y.(^)|-^/ 



^ 




Ml 


87r2 ^ ' 




Ml 


= ^MM^ 
87r2 ^ ' 




Ml 



7. + 1 + - ^-§T! (74) 



c 
,2 



where fc' is defined as k' = (ik and we have used the result 



°° k' 



cifc'-;^ — = — (75) 
+ 1 12 ^ ^ 

Thus, the critical temperature for the chiral symmetry restoration is given by 



which shows that the critical temperature is characterized by the quadratic term evaluated 
in the LORE method, which differs from the dimensional regularization where the quadratic 
term is in general suppressed. 

So far, we have explicitly shown the mechanism of dynamical spontaneous chiral symme- 
try breaking and its restoration at finite temperature. 

We are now going to present the explicit expressions for the masses of the scalar mesons, 
pseudoscalar mesons and/or light quarks. To be manifest, let us first write down the scalar 
and pseudoscalar meson matrices 

V2^-(*t* (77) 

and 

v^n= ^0 , (78) 

V ^ "71 + 73/ 
Keeping to the leading order of the current quark masses, we have 

ml.AT) = ml{T) ^ ^^(-« + = (79) 

for the pseudoscalar mesons, and 

m\.^{T) = ml{T) ^ 3(m^(T) + ml{T)) = Qm\T) (80) 

"0 

for the scalar mesons. Where /Xp is given by 

^^l{T) = mliT) + 2\{T)vl{T))vo{T) = ^i]{T)vo{T) = ^,vl{T) (81) 

where we have used the minimal condition Eq. fl66p and the relation Eq. fl73p . 

Note that in obtaining the above results for the scalar and pseudoscalar meson masses, 
the SU(2) triplet and singlet mesons have the common masses: m^^ = and = m^/, 
which is the refiection of the present assumption of the exact U{1)a symmetry. However, 
as we discussed in Sec. 2, in the real world such a symmetry gets quantum anomalous and 
will be broken down by the instanton effects, which is ignored in our present consideration. 



V. PREDICTIONS WITH INPUT PARAMETERS AT LOW ENERGY AND 
CRITICAL TEMPERATURE IN CTDM 



In order to make numerical predictions for the temperature dependence of the masses of 
the hght scalar and pseudoscalar mesons, it needs to fix the values of input parameters in 
the effective chiral Lagrangian with finite temperature. There are in general five parameters: 
/ij (/^m)) l^li "^o, and a universal current quark mass m. To fix the parameters, we shall 
use the constraints at low energy with zero temperature. 

In general, the minimal condition Eq.f l6ip with different quark masses will lead to two 
constraints by expanding the equation with respect to the current quark mass up to the 
order of m? . For the equal mass case, we get the following minimal condition 

Jx) + 2\vl = (82) 

with 

Tl'f-W-4)-^Ml-r)+.|^l (83) 

Note that in obtaining the above result one needs to keep to the order of in the current 
quark mass expansion. 

As we have shown in Sec. II that in order to have well-defined QCD current quark masses, 
it requires that 



^-1)M = M, i.e. /i^ = 2;uj (86) 



which fixes the parameter 



/So = 4 = 2 (87) 

Also, from the original Lagrangian of chiral dynamical model Eq.(|3]), the auxiliary fields 
$jj are found to be given by the quark fields as follows 

= jWLi + ^Mi, (88) 

By assuming that the quark condensation is almost flavor independent, i.e., < uu >~< dd >, 
and combining the condition Eq. (186!) , then the dynamical quark masses take the simple form 

rh = v — m = Vo + {Po — 1)"^ = Vo + m (89) 

which may be identified with the expected constituent quark masses after dynamically spon- 
taneous symmetry breaking, and Vo is caused by the quark condensation 

Vo = -7^{qq), q = u,d (90) 



To determine the remaining parameters, we consider the following constraints. There are 
two constraints arising from the pseudoscalar sector. One is from the normalization of the 
kinetic terms. 

A^o' = f (91) 
After some manipulation, the equation can be transformed into the following form 

^0(^)^0 = - A? - {UmeVf (92) 



where we have used the pion decay constant /tt — 94MeV. 

The other comes from the current quark mass and pion mass via the relation Eq.( !79|) . 

20 
20 

+ 0.08^ 



Taking the pion mass m^o,± ~ 139 MeV and the VEV Vo = 350+2° MeV or alternatively the 
current quark mass m = 4.76To ogMeV as the inputs, we obtain the following relation 



2 f2 

2^i}vo = -m = ^^^^^ = {2mtlMeVf (93) 



With the above relations and constraints Eqs.([82]), ([83]), ([86]), ([90]), ([92]), ([93]), all the 
parameters can be completely determined 

Vo ~ 350+^OMeK 

32 /\/r^T/ ,, Ql 0-3 



Mc ~ 88l;^^MeV, /i, ~ 312;^iMe1/ 
4 



= 2/iJ = (226;?Me\/)2 



/3o = 2, 7o = ^ = 0.209;°:°!} 



V. 



(qq) = ~{262tlMeVy (94) 

With these parameters, we can immediately obtain the critical temperature for the chiral 
symmetry restoration 



Tc = \l^2\^c - /^^In ^ - 7. + 1 + l/2(^))] ^ 2mil,MeV (95) 



which is consistent with NJL model prediction |20H23 . 



VI. CHIRAL SYMMETRY RESTORATION AND CRITICAL PHASE TRANSI- 
TION OF LOW ENERGY QCD 

In this section, we will present numerical predictions based on the CTDM. Especially 
we will show the thermodynamic behavior of the VEV fo(T), pion decay constant fjriT), 
the quark condensate {qq){T) and the masses of pseudoscalar mesons m^o,±(T). Their 
temperature dependence and the properties of critical phase transition are plotted in all the 
diagrams with adopting the central values of the quantities listed in Eq. (lM]) . 

From the gap equation Eq.( l72p . we can numerically solve the vacuum expectation value 
Vo{T) at any finite temperature until the critical temperature where Vo{T) approaches to 
vanish. The result is shown in Fig.([T]). 
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FIG. 1: Temperature dependence of VEV Vo 



By the normalization of kinetic terms of pseudoscalar sector Eq. fl91l) . we can obtain the 
expression determining the pion decay constant at finite temperature 

MT) = ^4X{TKiT) = 2^o(r)^^Zo(r) (96) 
which is presented in Fig. ([2]) 



fAT) 




FIG. 2: Temperature dependence of the pion decay constant 



Furthermore, the quark condensate {qq){T) is given by: 

(gg)(T) = -2fijiT)voiT) = -2^ovl{T) (97) 

its variation with respect to temperature is displayed in Fig. ([3]) 

The leading order approximation of the pseudoscalar meson mass m^o,± with respect to 
current quark mass m are expressed in Eq. fl79l) . which is shown in Fig.(jl]). 
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FIG. 3: Temperature dependence of the quark condensate 
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FIG. 4: Temperature dependence of the pion mass 



Let us now turn to thermodynamic property of the pseudoscalar meson mass Eq.( ff9l) : 



4/i3,(T) AvoiT)fx}iT) 



PAT) 



-m 



PAT) 



-m 



X{T)vo{T) 



m 



(98) 



which exphcitly shows that when keeping the mass scale /i^ to be a temperature-independent 
constant, the thermodynamic mass of the pseudoscalar meson becomes divergent near the 
critical temperature Tc as Vo(Tc) = 0, which is obviously contrary to our intuition. This 
is a manifest reason why we should make an assumption for the temperature dependence 
of /ij(T) given in Eq.f l73p . which can lead to the expected thermodynamic behavior for the 
pseudoscalar meson mass near the critical point 



(99) 



which is shown in Fig.(jl]). 

According to the above derivation, we see that the phase transition of chiral symmetry 
restoration is second order in our simple modeL Thus, it is natural to determine the critical 
behavior of all the quantities discussed previously. Now we would like to find the the scaling 
behavior of the vacuum expectation value Vo{T) near the critical point. By expanding our 
gap equation Eq. fl72l) around the critical temperature Tc with respect to the small value of 
VEV Vo{TY up to the order of fo(T)^, we can obtain: 

\n,{T^^ - T^) - CvoiTf = 0, (100) 
o 

where C = ^7(0, + | — 7o and 7(s,x) = f^'^e'^dt is the lower incomplete gamma 
function. Given above equation, we can easily obtain: 

^o=y^(Ti^-T^)^oc(T.-T)i (101) 

Thus, the critical dimension is /3 = 0.5. Other quantities such as f-n{T)^ m^o,±(T) and 
{-{qq){T)f'^ all have the same scaling behavior. Such a critical behavior is not accidental, 
which can actually be understood from the fact that near the critical temperature all these 
quantities are proportional to the VEV Vo{T) with A(T) keeping fixed to A(Tc) 7^ 0. 



VII. CONCLUSIONS AND REMARKS 

In this paper, we have extended the chiral dynamical model to the chiral thermodynamic 
model by adopting the CTPGF approach. The resulting effective chiral Lagrangian for 
the composite meson fields at finite temperature is similar to the one of the CDM, but 
all the couplings and mass scales become temperature dependent. We have discussed in 
detail the finite temperature behavior of CTDM. Much attention has been paid to the 
thermodynamic chiral symmetry breaking and its restoration at finite temperature. After 
fixing the free parameters in the effective chiral Lagrangian, we have determined the critical 
temperature for the chiral symmetry restoration, its value has been found to be around 



Tr ~ 200 MeV which is consistent with other predictions based on the NJL model|20- 



23| . We have also explicitly presented the thermodynamic behavior of several interesting 
quantities which include the vacuum expectation value VEV Vo{T), the pion decay constant 
fniT), the quark condensate {qq)(T) and the pseudoscalar meson mass m^o,±(T), they all 
display the property of the chiral symmetry restoration at the critical temperature T^. From 
the numerical calculations, we have shown that they all have the same scaling behavior near 
the critical point. It is interesting to note that the mass scale Hf for the four quark interaction 
in the NJL model should be temperature dependent at finite temperature as expected from 
the gluon thermodynamics, its thermodynamic behavior near the critical point is required 
to be same as the one of the chiral symmetry breaking. As a consequence, we are led to the 
assumption that /ij(T) = 7of^(T) in order to yield the expected thermodynamic behavior 
of the pion meson mass and to avoid the divergent behavior near the critical point of phase 
transition. Finally, we would like to remark that as limited from our main purpose in the 
present paper we have only considered two flavor quarks and ignored the important instanton 
effects and U(1)a anomalous effect, which prevents us to discuss some other interesting 
properties, such as the large strange quark mass effects and the anomalous U{1)a symmetry 
restoration at finite temperature, we shall investigate those interesting effects elsewhere. 
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Appendix A: Brief Outline on Closed- Time- Path Green Function (CTPGF) For- 
malism 



The formalism used in zero-temperature quantum field theory is suitable to describe 
observables (e.g. cross-sections) measured in empty space-time, as particle interactions in an 
accelerator. However, at high temperature, the environment has a non-negligible density of 
matter which makes the assumption of zero-temperature field theories inapplicable. Namely, 
under those circumstances, the methods of zero-temperature field theories are not sufficient 
any more and should be replaced by others, which is closer to thermodynamics where the 
background state is a thermal bath. Therefore we shall develop quantum field theory with 
finite temperature which is extremely useful to study all phenomena due to the collective 
effects, such as: phase transitions, early universes, etc. There are several approaches for 
the finite temperature field theories, in this appendix we shall focus on the closed-time-path 
Green function (CTPGF) formalism which is simply applicable in our case. The CTPGF 
formalism, developed by Schwinger [13] and Keldysh |l5|, has been used to solve lots of 
interesting problems in statistical physics and condensed matter theory It is generally 
believed that this technique is quite efficient in investigating the nonequilibrium and finite 
temperature dynamical systems, this is because such a formalism naturally incorporates 
both the statistical and dynamical information 16, 17|. Excellent review articles [ItI . 2?! 
and monographs [isl, 19[ have described different aspects of these issues. In this appendix, 
we will briefly outline the main method with the Schwinger-Keldysh propagators and the 
universal Feynman rules for the general theory. 

Let us begin by the general discussion of statistical physics. A dynamical system can be 
characterized by its Hamiltonian H and a statistical ensemble of this system in equilibrium 
at a finite temperature T = (in units of Boltzmann constant) is described in terms of a 
partition function 

Z{(5) = Tr p(/3) = Tr e'^^ (Al) 

Here p(/3) is known as the density matrix operator and H can be thought of as the generalized 
Hamiltonian of the system. For example, for a canonical ensemble in which the system can 
only exchange energy with the heat bath, is defined as: 

H = H (A2) 

while for a grand canonical ensemble in which the system can not only exchange energy with 
the heat bath but also exchange particles with the reservoir, "H is taken as 

H = H-i2N (A3) 

where /i is the chemical potential and N represents the particle number operator. 
A observable in a statistical ensemble is the ensemble average for any operator 

{0), = ^^Ttp{P)0 (A4) 



Since the partition function and ensemble averages involve a trace operation, this feature 



leads to the famous KMS (Kubo-Martin-Schwinger) relation. 



1 

m 
1 



Tre-^^Ci(t)02(t') 



(A5) 



Note that KMS relation only rely on the trace operation and does not depend on any 
periodicity property of operators along the temperature (imaginary time) interval. 

Now it is easily seen that the operator e~'^^ in the definition of the partition function is 
very similar to the time evolution operator in the imaginary time axis e~*(~*^)^ j29|. We 
promote this similarity and analytically extend the time variable to the complex plane. So 
the operator e~'^^ would live on the line interval which is parallel to the negative imaginary 
time-axis, with its length /3. The analogy implies that we can define our theory on some 
certain contour on the complex t-plane. 

The contour should satisfy the following several requirement: (i) The two endpoints of 
the contour must be fixed in an interval in a line parallel to the imaginary axis with its length 
/3. The stating point A (corresponding to time ti) and the ending point B (corresponding 
to tf = U — i(3) are identified, and one requires that 0\b = 0\a if O is bosonic and and 
0\b = —0\a if O is fermionic; (ii) For a system whose spectrum of the Hamiltonian is semi- 
positive (at least bounded below due to the stability of the system), the contour needs to 
have a monotonically decreasing or constant imaginary part for the reason of the convergence 
of the complete partition function. The same result can also be obtained by analyzing the 
convergence of the two-point Green function. 2^; (iii) The contour needs to pass the whole 
real axis of t-plane on which the field operators Oi{t) are defined (Real-time Formalism. 
Otherwise, like the imaginary time formalism, the field operators need to be analytically 
extended to the imaginary axis first). 

The particular family of such real time contours is depicted in Fig. |5] where the contour 
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FIG. 5: Contour used in the real time formalism 



C is C = Ci U C2 U U C4. The contour Ci goes from the initial time ti to the final time 



tf, C3 from tf to tf — ia, with < a < /3, C2 from tf — ia to ti — ia, and C4 from ti — io" to 
ti — ip. Different choices of a lead to an equivalence class of quantum field theories at finite 
temperature. Our preferred choice is the Schwinger-Keldysh one with a = 0. 

With this specific contour, the action of a field configuration is the sum of contributions 
from the three parts, 

tf t{ 

S = jdtL{t) = j dt L{t) - j dt L{t) - i J dr L{ti - ir) (A6) 
c ti ti 

where 

L{t) = j dxC[(t){t,x)] , (A7) 

and C is the Lagrangian density. In the following we will take the theory of a scalar field 
as an example. However, in the limit ti — t- —00 and tf ^ 00, it can be shown that the third 
branch gets decoupled from the other two (the factors in the propagators connecting such 
branches are asymptotically damped). Consequently, in this limit, we are effectively dealing 
with two branches leading to the name " closed time path formalism" . In this contour, then, 
the time integration has to be thought of as 



00 



dt= I dt+- dt^ (A8) 



00 



where the relative negative sign arises because time is decreasing in the second branch of 
the time contour. 

The advantage of introducing the contour C is that one can introduce the sources coupled 
to the field which is not vanishing on the two Minkowski parts of the contour. This 
procedure would give us the generating functional 

(00 00 \ 

iS + i J dt^ J dx Ji{x)(j)i{x) ~^ j J dx J2{x)4)2{x) | . (A9) 
—00 —00 / 

Here Ji^2 and 0i^2 are the sources and fields on the two Minkowski parts of the contour, i.e., 

Ji{t,x) = J(t^,x) , 0i(t, x) = 0(t+, x) , (AlOa) 

J2(t, x) = J{t_,x) , 02(i, x) = <p{t_,x) . (AlOb) 

By taking second variations of Z with respect to the source one finds the Schwinger- 
Keldysh propagator 

1 5HnZ[Jr,J2] .f Gu -Gu\ 
^^'^^(^ -y^ = -^ 5Ux)5My) = ' [ -G2, G22 ) ■ 

In the operator formalism, the Schwinger-Keldysh propagator corresponds to the contour- 
ordered correlation function. In the single time representation this means: 



iGu{t,x) = (T0i(t,f)0i(O))/3 , iGi2(t,x) = (02(O)0i(t, x));? , 

zG'2l(t, X) = {(t>2(t, x)0i(O))/3 , iG22(t, X) = (T02(t, x)02(O))/3 . 



where T (T) denotes normal (reversed) time ordering, and 



= e^^*-*^->(0)e-*-^*+^^-^, (A13a) 
(f)2{t,x) = e^^(*-^'^)-^-P->(0)e-^^(*-^")+^-^-^\ (A13b) 

Let us now consider the free real scalar theory. If one goes to the momentum space and 
by inserting the complete set of states into the definitions flA12p . one finds the explicit form 
of the previously defined Schwinger-Keldysh propagator. 

^Giiik) = + 2nnB{uj)6{e - m^) , u = \ko\ (AM) 

iGu{k) = 27r[nB{co) + e{-ko)]6{P - m^) , (A15) 
iG2i{k) = 27i[nB{uj) + e{ko)]S{k^ -m^), (A16) 

iG22ik) = + 27mBico)5ik^ - m^) . (A17) 

f^Z _ ^Z _j_ 



Or in the matrix form: 

^Gp{k) 



2n9iko)5ik' - m') -p^fc; 
+ 27cnB{co)6{e - m^) J M (A18) 



where nB{uj) = stands for the Bose-Einstein distribution function. Note that the 

propagator is a 2 x 2 matrix, a consequence of the doubling of the degrees of freedom. 
However, the propagators (12), (21) and (22) are unphysical since at least one of their time 
arguments is on the negative branch. They are required for the consistency of the theory. 
The only physical propagator is the (11) component shown in Eq. flA14p . 

For perturbative calculations we need to know the complete Feynman rules besides of 
the propagators. From the generating functional Eq. flAQI) and the action Eq. flAGD defined on 
contour C Eq. flASp . we see that the complete theory contains two types of vertices- type-1 for 
the original fields 0i(x) while type-2 for the doubled fields 02 (2;)- The vertices for the partner 
fields will have a relative negative sign corresponding to the original vertices, because time is 
decreasing in the negative branch. The four possible propagators, (11), (12), (21) and (22) 
defined above connect them. All of them have to be considered for the consistency of the 
theory. The golden rule is that: physical legs must always be attached to type 1 vertices [28l]. 
For other Feynman rules, including the integration measure, the symmetry factors involved 
in Feynman diagrams, the topology of the Feynman diagrams, etc. are all the same as the 
zero-temperature field theory. 

For the application to the present paper, we also need to know the Schwinger-Keldysh 
propagator for fermions as the quark fields here are represented as the chiral fermion fields 

27r^(-A;o)5(P - 



2 



m 

2 _2 ^ 



J fc2- 

1 r 



-27rnp{u)6{k' -m')\^ ^ ^ j] (A19) 
where npiuj) = -. stands for the Fermi-Dirac distribution function. 



When transforming into the Euchdean spacetime the Schwinger-Keldysh propagator de- 
fined above becomes 

-2^inF{uj)5{kl + m^) 1^ J | ^ ] (A20) 

In Sec. 3 the factor —i is canceled by the factor i in the integration measure transformation 
d^k — )■ id^kE- 



Appendix B: Derivation of Chiral Dynamical Model in the Chiral Rotated Basis 

In this appendix, we will derive the effective chiral Lagrangian for mesons in the so- 
called chiral "Rotated Basis" 13|]. Although the obtained effective chiral Lagrangian will not 
change, it is more transparent to see the chiral symmetries and their spontaneous breaking 
in this derivation. Let us begin with the effective Lagrangian 

^^IfMi q) = qi^id^q + qil^^A^qL + qnl^A^i^qR - [qd^ - M)qR + h.c] (Bl) 
where the auxiliary meson fields is defined as in Eq.([S]) 

J- J- n .2n(3:) 

$(x) = U{x) ^ = ai^) = 

3 3 

(P\x) = (f){x) = J2 <P\x)T'', ^\x) = n(x) = J2 ^''(x)T^ (B2) 

a=0 a=0 

where 11 (x) and (f){x) represent the pseudoscalar and scalar mesons respectively. Note that 
except for the mass term or the source term, the Lagrangian is invariant under the trans- 
formation of the local U(2)l x U{2)fi chiral symmetry: 

qM -> gLix)qLix), qnix) gR{x)qR{x)] ^{x) gL{x)^{x)g'^^{x), 

^L/. -> g\ALt,gL{,x) - igldf^giix), Ar^ ^ 9RAR^,gR{x) - igRd^gRix), (B3) 

The transformation for can also be written in terms of the fields and as: 

^ h{x)(t){x)h^{x), ^l(x) -> gLix)^Lix)h^{x) = /i(x)^l(x)c/]j(x). (B4) 

Let us now introduce new quark fields, which is referred to the chiral "Rotated Basis" in[l3|. 

qL = ^lQl, qL = Ql^I, 

qR = ^IQri qR = Qr^l- (B5) 

With this new quark basis, we can rewrite the Lagrangian Eq.f lBip in the following form 

jC'^ffiQ, Q) = QYtd^Q + QlI^L^Ql + Qr^R^^Qr - [Ql{'P - M)Qr + h.c], (B6) 
where the fields L^, and Ai are defined as 

= eML^L + ^ei<9^a, R^ = ^lAr^^I + t^Ld^^l 
M ^ ClMCi, ^ iLM^iL. (B7) 



In the above "rotated basis", the quark fields Ql{r){x) transform only under the diagonal 
f/y(2) symmetry: 



Ql{x) ^ h{x)QL{x), Qr{x) -> h{x)QR{x). (B8) 

Thus, the quark fields Qjar){x) are much like the "constituent quark" defined in the non- 
relativistic quark model When the chiral symmetry is spontaneous breaking and the 
meson field acquires the vacuum expectation value (VEV) < >= V, Ql{r) "will 
obtain a mass term QiiY — M.)Qr + h.c. In the case where each quark fiavor possesses 
the universal current mass m, the VEV matrix is diagonal V = v ■ I, here / is the identity 
matrix in the fiavor space. The mass term is {v — m)Q iQ r + h .c. , namely the mass of quarks 
Q{x)l(r) is the dynamical quark mass fh = v — m defined in SecJITl 

Note that the Lagrangian Eq. (lB6p has the same structure as the original one Eq. (lBip 
except for the definition of the mass and gauge fields. Thus, we may expect that the 
effective chiral Lagrangian for the meson fields has the same structure as Eq. fHSj) . Indeed, 
by integrating over the quark fields following the procedure from Eq.Q to Eq.f pOj) . we 
obtain 

S^Re =YJ ^'^^ J ^^^''^ 1^^^' + - 1^ ^^9) 

where A'^ above is defined as 

A'e = U^Pr + 0^0Pl - 11 ■ D'e<PPl - tl ■ D'E(j)^Pn 

-^,^^'r,uPl - ^.u^^uPr + {^D'^,){^D'^,) + 2 A; ■ {iD'^), (BIO) 

and 



= id^cj) + - (BU) 
iD'e^ = id^ + R^Pl + L^Pr, (B12) 

4) = (t)-M. (B13) 

In order to derive the effective action for meson field, we redefine A^' = A;^ + A'^ to the 
following two terms: 

A^' = A;2 + A'^ = Ao + A'^, (B14) 

with 

Ao = fc2 + M\ 

A'e = [(00^ - MM^)Pr + (0^0 - M^M)Pl - 11 ■ D'j,</)Pl - 11 ■ D'j,<P^Pr 
-^,.J"r,.Pl - ^,uJ"l,uPr + i^D'E^){^D'^^) + 2k ■ (tD'^)] 
+{MM^ - M^)Pr + {M^M - M^)Pl, (B15) 

where M is supposed vacuum expectation values (VEVs) of $, i.e., < $ >= M which is real 
and Ai = .^]^M^]^. As we will see that the terms in the third line of the definition of A'^ is 
crucial to prove the equivalence of the obtained effective action for mesons in this rotated 
basis to the one given in Eq.f HSj) . 



Now if we regard A'^; as the perturbation and take Zq = (det Aq) 2 as before, we can 
expand the effective action Eq. flB9|) according to A'^: 



If we only keep the leading two terms in the expansion as in the previous sections, we can 
obtain: 

2 J ""J (27r)4 ^^L^o 2Al 



= 

167r2 

-^Lo[D^0 • D'e^^ + D'e^^ ■ D'e^ + _ MM^f + (0^0 - M'^Mf] 
+MlL2\{MM^ - M^) + {M^M - M^)] 

-iLo[(00^ - MM^){MM^ - M^) + (0^0 - M^M){M^M - M^) 

+{MM^ - M2)(00t _ MM^) + {M^M - M^){^^ - M^M) 

{MM^ - M^f + {M^M- M^f]}, (B17) 

where the matrix Lq and L2 are defined in Eq. flH^ . It is easy to see that the last four lines 
of terms vanish if the different flavors of quarks have the same current mass which leads to 
the same vacuum expectation value for each flavor: 

MM^-M^ = (eiMei)(aMeL)-M' 

= ilM'^iL - = 0. (B18) 

The last line is valid since the mass matrix M is diagonal with the same eigenvalues and M 
can commute with the SU(2) matrix ^l. 

Note also that in the case of 2 flavors with the universal current quark mass, after the 
chiral symmetry SUl{2) x SUr(2) is broken to SUyi^) and the field 0(a;) acquires vacuum 
expectation value (VEV) = f ■ /, we have: 

= {v- m)ii{{DE^JJ)i{ = -{v - m)iiL{DE,U^)iL. (B19) 

In order to prove the last two equalities, we have to use the definition of f/ = and the 
identity d^j^^i = —C,L{dfj.^l)^L- Therefore, from the kinetic terms for 0(x) in Eq. (]B17p . we 
can obtain the kinetic term for the pseudoscalar meson field U = e'^'^(^)/-^'" : 



Nc(v — m) 



J d^XEtrE[Lo{DE^U){DE^.U^)]. (B20) 



However, when the chiral symmetry is restored, as discussed in the context of CDTM, this 
term will disappear as the VEV of 0(x) vanishes. Thus, the discussion in this "rotated 
basis" gives a more transparent picture of Goldstone boson character of the pseudoscalar 
mesons. 

Next we would like to prove the equivalence between the effective chiral Lagrangians 
Eq. (H3l) and Eq. (lB17l) . According to the definition of the transformations Eq. (1B7l) . we can 
easily obtain the following relations: 

D'e,^ = i{{DE,^)i[ D'e^4>^ = UDE,^^)iL, (B21) 

where De^i^ and DEfi^"^ are defined as in Eq.( |22l) . Thus, the first two lines of terms in 
Eq. (lB17p can be written in a form with respect to $: 



SMne = yI d'xEtTE{M'M^[iU^-M')^, + U^^^-M')^l] 
-hoi^liDE^ ■ De^^)U + UDe^^ ■ De^)^1 

= yI d'xEtTE{M'MLL2^l)i^^^-M') + i^iL,^L)i^^^-M^^^ 

+ {^[Lo^l)[De^^ ■ De^ + i^^^ - M')']}}. (B22) 

In general the matrices C,l and Lq (L2) do not commute with each other when different flavors 
do not have the same current masses. So the part of effective chiral Lagrangian shown in 
Eq. flB22p is in general not equivalent to Eq. fH3|) by just comparing to the same truncated 



terms. In fact, by taking into account the higher-order terms in Eq.( ]B16l) and the last four 



lines of terms in Eq.( ]B17|) . it is expected that the extra terms would cancel the unwanted 
terms in Eq. flB22p due to non-commutativity of C,l and Lq {L2). In our present case, as we 
only consider two flavors with the same current quark mass, C,l can commute with Lq and 
L2, there is no such extra terms. Furthermore, as mentioned before, the last four lines of 
terms in Eq. (lB17p will vanish due to commutativity of C,l and M. Thus, we arrive at the 
following effective chiral Lagrangian: 



-^Lo[De^ ■ De^^ + De^^ ■ De^ + ($$"^ - M'^f + (l>"f<l> - M'^f]], (B23) 



which is exactly agree with Eq.f l43p obtained in the original unrotated basis. 

Following the procedure in Sec JIIIl we can derive the effective chiral Lagrangian at finite 



temperature in the chiral thermodynamic model (CTDM), which is: 
S^ER, - Yq^I d^XEtrF{M'^L2{T)[{<l)<l)^ - MM^) + (0^ - M^M)] 

-^Lo{T)[D'J ■ D'e^^ + D'e^^ ■ D'e^ + _ MM^f + (0^0 - Mf] 
+MlL2{T)\{MM^ - M^) + {M^M - M^)] 

-^Lo(T)[(#t _ MM^){MM^ - M^) + (0+^ _ M^M){M^M - M^) 

{MM^ -My + {M^M-My]}. (B24) 

By using the same argument as the one for the zero-temperature chiral dynamical model 
(CDM) analyzed in the text, we can prove the equivalence between Eq. flB24p and Eq.f lSU]) . 
Thus, based on the equivalent Lagrangians, the resulting physics, especially the spectrum 
of mesons, will not be changed. 
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